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1 Based on observations collected at the European Southern Observatory, La Silla, Chile, using SOFI at 
the 3.5m NTT, within the observing programs 64.N-0038 and 68.D-0287 
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We present deep J,H,Ks photometry and accurate Color Magnitude Diagrams 
down to K ps 18.5, for a sample of 13 globular clusters in the Large Magellanic 
Cloud. This data set combined with the previous sample of 6 clusters published 
by our group gives the opportunity to study the properties of giant stars in 
clusters with different ages (ranging from ~ 80 Myr up to 3.5 Gyr). Quantitative 
estimates of star population ratios (by number and luminosity) in the Asymptotic 
Giant Branch, the Red Giant Branch and the He-clump, have been obtained and 
compared with theoretical models in the framework of probing the so-called phase 
transitions. The AGB contribution to the total luminosity starts to be significant 
at ps 200 Myr and reaches its maximum at 5 — 600 Myr, when the RGB Phase 
Transition is starting. At ~ 900 Myr the full developing of an extended and well 
populated RGB has been completed. Both the occurrence of the AGB and RGB 
Phase Transitions are sharp events, lasting a few hundreds Myr only. These 
empirical results agree very well with the theoretical predictions of simple stellar 
population models based on canonical tracks and the fuel-consumption approach. 

Subject headings: galaxies:star clusters - Magellanic Clouds - infrared: stars - 
techniques: photometry 



1. Introduction 

The red giant star populations play an important role in the spectral evolution of a 
simple stellar population (SSP), dominating the bolometric luminosity of a SSP after a few 
hundred million years. Two special events (the so-called phase transitions, hereafter Ph-T) 
due to the sudden appearance of stars with C-0 and He degenerate cores mark the spectral 
evolution of a SSP. The first event corresponds to the appearance of bright Asymptotic 
Giant Branch (AGB) stars; the second correspond to the full development of the Red Giant 
Branch (RGB). Several theoretical models (Renzini & Buzzoni 1986; Bruzual & Chariot 
1993; Maraston 1998) investigated the impact of the AGB and RGB stars on the total light 
of a SSP. The AGB Ph-T is predicted to occur after ps 10 8 years and the RGB Ph-T after 
ps 6 • 10 8 years. These events must be empirically calibrated by using suitable templates 
of SSPs, in order to use them for predicting the age of more complex stellar systems like 
galaxies. The Large Magellanic Cloud (LMC) globular cluster system represents a gold-mine 
to study the AGB and RGB properties with varying the age and chemical composition of 
the stellar population. Indeed, these clusters cover a wide range in age (Searle, Wilkinson, 
& Bagnuolo 1980; Elson & Fall 1985; Girardi et al. 1995; Brocato et al. 2001), in metallicity 
(Sagar & Pandey 1989; Olszewski et al. 1991; Hill et al. 2000) and in integrated colors (van 
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den Bergh 1981; Persson et al. 1983), providing the ideal template to investigate the spectral 
behavior of a SSP. 

The AGB and RGB phases are intrinsically difficult to sample with the necessary statis- 
tical significance, because of their relatively short lifetimes. Several searches investigated the 
ratio between C and M stars and its correlation with metallicity in the LMC field (Blanco 
et al. 1980; Cioni & Habing 2003; Cioni et al. 2005), finding more C stars in lower metal- 
licity environments. Frogel et al. (1990, hereafter FMB90) presented a detailed study of C 
and M-type AGB and their contribution to the cluster total light as a function of Searle, 
Wilkinson, & Bagnuolo (1980, hereafter SWB) classes. FMB90 assumed that the separation 
between AGB and RGB is positioned at M{, ;=-3.6 and stars brighter than that value are 
on the AGB. They found C-stars only in the intermediate age SWB IV- VI clusters. The 
SWB VII clusters are too old, and their AGB stars have too low mass stellar progenitors 
for evolving as C-stars, while SWB I-III clusters are still too young to have stars evolving 
along the AGB. In these objects the brightest red stars are M supergiants. For intermediate 
age SWB V-VI clusters 40% of the bolometric luminosity is contributed by AGB stars. In 
particular C-stars, when present, account for 50-100 % of the bolometric luminosity of the 
AGB. 

The observational results by FMB90 have been used by Maraston (1998, 2005) to calibrate 
the energetics and compositon of the TP-AGB phase to be inserted in stellar population mod- 
els. It is therefore very interesting to check those observational results with new, accurate, 
data. 

The first observational evidence of the RGB Ph-T in stellar clusters has been presented 
by Ferraro et al. (1995, hereafter F95) who analyzed the red giant contents of 12 intermediate 
age LMC cluster. This survey was performed with the first generation of IR imagers, and 
sampled only the brightest portion (K<14.3) of the RGB. Despite of its limited luminosity 
sampling, this pioneering work demonstrated that near-IR observations are indeed crucial 
in order to study the photometric properties of AGB and RGB sequences, providing the 
highest sensitivity to the physical parameters of cool stars. 

Recently, the superior performances of the the 1024x1024 IR array detector of SOFI at 
the ESO/NTT, allowed to measure the entire extension of the RGB for a complete sample of 
LMC clusters spanning a wide range of ages. Ferraro et al. (2004, hereafter F04) presented 
deep J,n,Ks photometry and accurate Color Magnitude Diagrams (CMD) down toKs; 18.5, 
i.e. ~ 1.5 mag below the red He-clump, for six intermediate age clusters (namely NGC 1987, 
NGC 2108, NGC 2190, NGC 2209, NGC 2231, NGC 2249). A quantitative estimate of the 
population ratios (by number and luminosity) between Red Giant Branch and He-clump 
stars for each target cluster, and comparison with theoretical predictions suggest that the 
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full development of the RGB should occur at t «j 700 Myr and be a relatively fast event 
(«5t « 300 Myr). 

The present paper shows the results for 13 additional LMC clusters, extending the 
discussion about the RGB Ph-T presented in F04 and studyng the contribution to the total 
luminosity of the AGB. The paper is organized as follows: Sect. 2 describes the observations 
and photometric analysis, Sect. 3 the near-IR CMDs. Sect. 4 present the basic assumptions 
on reddening and the adopted age calibration. Sect. 5 discusses the procedure adopted to 
derive complete and accurate star counts and luminosity contribution. Sects. 6 and 7 discuss 
the phase transitions, while in Sect. 8 we draw our conclusions. 



2. Observations and photometric analysis 

J,H,Ks images of 13 globular clusters in the LMC (see Table 1 have been obtained 
at ESO, La Silla, on December 28-30, 2001, with the near IR imager/spectrometer SOFI 
(Moorwood, Cuby & Lidman 1998) mounted at the ESO 3.5m NTT. SOFI is equipped with 
a 1024 x 1024 Rockwell IR- array detector. All the observations presented here, have been 
performed with a scale of 0.292"/pixel, providing a ~ 5' x 5' field of view, each frame. 
The observations were obtained in good seeing conditions (0.7" — 0.8" on average). Total 
integration times of 2 min in J, 4 min in H and 8 min in Ks split into sets of shorter exposures 
have been secured, allowing to obtain accurate (S/N>30) photometry down to J~19 and 
H,Ks^l8.5. 



l : J 



Clusters have been selected accordingly to the s-parameter, defined by Elson & Fall 
(1985) as a curvilinear coordinate running along the mean locus defined by the (U-B) and 
(B-V) integrated colors of the LMC clusters, and directly linked to the cluster age. The 
6 clusters presented by F04 have s— parameter ranging between 34 and 37. The database 
presented here include 7 additional clusters with s = 23 — 31 and 6 with s = 37 — 45. Table 1 
lists the basic properties of the clusters discussed in the paper: s values (Elson & Fall 1988) , 
E(B — V) (Persson et al. 1983) and metallicity (Sagar & Pandey 1989; Olszewski et al. 1991; 
Hill et al. 2000; Dirsch et al. 2000). Indicative [Fe/H] values vary by a factor of ~10, between 
-1.2 and solar with an average value of ~ —0.6 ± 0.1, but precise abundance determinations 
via high resolution spectroscopy still lack for the majority of the LMC clusters. 

A control field a few arcmin away from each cluster center has been also observed using 
the same instrumental configuration, in order to construct median-averaged sky frames and 
for decontamination purposes. A large sample of high S/N flat-fields in each filter have been 
acquired by using a halogen lamp, alternatively switched on and off. The final cluster and 
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control field frames have been sky-subtracted and flat-field corrected. For each cluster, the 
dataset includes three J,H,Ks images, centered on the cluster and three corresponding frames 
of an adjacent field. 

The photometric analysis was performed by using daophot-II (Stetson 1987). For each 
observed field all the images in the J,H,Ks filters were carefully aligned and trimmed in 
order to have three output images, one per filter, slightly smaller than the original ones 
but perfectly registered. Then, a daophot-II Point Spread Function (PSF) fitting run was 
applied to each image. The output catalog with the instrumental magnitudes has been 
checked for any spurious detection or missing object (typically 3-4 stars at most) which 
have been included in the catalog by hand. 

Each catalog typically contains 1000-1200 stars, with the exception of 3 poorly populated 
clusters counting less than 700 stars (NGC 1651, 2162 and 2173) and 2 rich clusters with more 
than 1900 stars (NGC 1783, 1978). The instrumental magnitudes have been transformed into 
the Two-Micron All-Sky Survey (2MASS) photometric system, by using the large number of 
stars (typically a few hundred) in common. The overall dispersion of these transformations 
is a < 0.01 mag in all the three filters. For sake of homogeneity, the output catalogs of F04 
(calibrated on the Persson et al. (1998) photometric system) have been also transformed 
into the 2MASS one, although the difference between the two calibration turns out to be 
negligible (~ 0.01 mag). 

The calibrated photometric catalogs in each filter were finally matched and merged 
together in a global catalog, using the CataXcorr and Catacomb softwares developed at the 
Bologna Observatory for an optimized cross-correlation. 



3. Color— Magnitude Diagrams 

Figs. 1-4 show the CMDs for the observed clusters and adjacent fields in the (K, (J— K)) 
plane. We have divided the clusters into two groups, accordingly to their s-parameter and 
different CMD morphology. The first group (see Fig. 1) includes clusters with s — 23 — 31; 
the second group (see Fig. 3) those with s = 37 — 45. Figs. 2 and 4 show the CMDs of the 
corresponding adjacent fields. x 

The CMDs of the first group of clusters (with s < 31) appear quite complex, hence 
particular care has been devoted to separate the cluster population from the LMC field. In 
order to help the reader identify the two populations we have plotted in the last panel a sketch 



1 It is evident the extreme field contamination for NGC 203f , located in the LMC Bar. 
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showing the mean location of cluster (grey regions) and field (dashed region) population in 
the CMD. The main properties of the CMDs shown in Fig. 1 can be summarized as follows: 

(i) Magnitudes as faint as K^ 18.5 have been measured. 

(ii) The brightest objects at K < 13 are likely AGB stars. 

(Hi) A blue sequence is clearly visible at -0.3< (J — K) < 0.3 and K>15.5, corresponding 
to the brightest end of the cluster MS. 

(iv) Helium clump cluster stars define a sequence at K = 13 — 14 and (J — K) ~ 0.2 
in the youngest objects (namely NGC2164, NGC2157, NGC2136). In older clusters (namely 
NGC2031, NGC1866, NGC2134, NGC1831), they define a clump at progressively lower 
magnitudes (K = 14 — 15) and redder colors (0.4 < (J — K) < 0.6). 

(v) The direct comparison of each panel in Figs. 1 and 2 clearly shows the significant 
contribution of the LMC field stars. As schematically shown in the last panel of Fig.l, most 
of the stars in the region 0.4 < ( J — K) < 1 and K > 12 are indeed LMC field stars, with a 
well defined RGB at 12 < K < 16 and the He-clump at K?»17. 

The observed CMDs for the second group of cluster show the following characteristics: 

(i) Magnitude limits down to Km 18.5, i.e. about 1.5 mag below the He-clump, which 
is clearly visible as a clump of stars at K ps 17. 

(ii) A well-populated and extended RGB. 

(Hi) The brightest objects at K < 12 are likely AGB stars. 

(iv) Unlike the first group of clusters, in these objects cluster and field populations do 
overlap. 



4. Basic assumptions 

4.1. Reddening 

Correction for extinction is computed accordingly to the E(B — V) values reported in 
Table 1 and the Rieke & Lebofsky (1985) interstellar extinction law. The infrared dust maps 
by Schlegel, Finkbeiner, & Davis (1998) in the direction of the observed clusters provide very 
similar (on average within ±0.03 dex) E(B — V) corrections, with the exception of NGC 2031, 
for which the discrepancy is about 0.1 dex. However, the overall impact of such a discrepancy 
on the infrared magnitudes is always small (well within 0.1 dex), hence reddening correction 
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is not a critical issue in this context. 

Absolute and bolometric magnitudes have been obtained by adopting a distance modu- 
lus (m — M) = 18.5 (van den Bergh 1998; Alves 2004) and suitable bolometric corrections 
by using the (J — K)o color and the empirical calibrations by Montegriffo et al. (1998). In 
computing luminosities, we adopted Mq° 1 = 4.74 and Mq = 3.41 for the Sun. In the 
following all the derived luminosities are expressed in unit of 1O 4 L . 



4.2. Age calibration 

A suitable calibration of the LMC cluster age is still a major concern since homogeneous 
determinations based on the MS TO, for a significant number of clusters, are not available 
yet. Following F04, here we use the s— parameter (see Sect. 2). Being a pure empirical 
quantity, it needs to be calibrated with age. As discussed in F04, the most used calibrations 
by Elson & Fall (1988) based on canonical models and by Girardi et al. (1995) based on 
overshooting models provide, somewhat surprisingly, very similar ages (within 10-15 %). 
Although a new calibration of the s-parameter as a function of age (see F04) is urged, in 
the following, as done in F04, we adopted the most recent one by Girardi et al. (1995) and 
based on the models with overshooting by Bertelli et al. (1994): 

logt = 6.227 + 0.0733 ■ s. 

We assume a conservative error of <5s=±l, which translates into a m 20% age uncertainty. 

Table 1 lists the s value for the entire sample of 19 clusters (13 presented here and 6 
discussed in F04, respectively) and their age derived from the quoted calibration. As can be 
seen the entire sample cover a large range of ages from 80 Myr to 3.5 Gyr. 



5. Star counts and integrated luminosities 

A quantitative analysis of the AGB and RGB populations (by number and luminosity) 
is crucial to empirically calibrate the relative lifetimes and to quantitatively evaluate the 
impact of each evolutionary stage on the total luminosity of a SSP. In order to obtain 
reliable stellar counts and luminosities in each branch, we proceeded as follows: (1) stars 
in each evolutionary stage have been identified on the basis of suitable selection boxes as 
defined in the CMDs (as shown in Fig. 6 and Fig. 9); (2) each sample of stars has been 
corrected for incompleteness, following the standard artificial star technique (see below); (3) 
the contamination from foreground/background stars in each population has been evaluated 



and statistically subtracted to the observed samples. 

While the definition of the selection boxes (step (1)) for the AGB and RGB populations 
are described in Sects. 6 and 7, respectively, in the following we briefly discuss the procedure 
adopted to perform step (2) and step (3). 



5.1. Completeness and field decontamination 

The degree of completeness can be quantified by using the widely-used artificial star 
technique (see the discussion present in Mateo 1988). For each cluster we have derived 
the RGB fiducial line and then a population of artificial stars, having magnitudes, colors, 
and luminosity functions resembling the observed distributions was generated and added 
to the original images. Since crowding effects are more severe in the central regions, the 
frame area sampling the cluster has been divided in three concentric regions (Region A, B 
and C, see left panel of Fig. 5) and the completeness has been estimated independently in 
each of them. The maximum spatial extension of each cluster has been estimated from the 
cluster radial density profile. A total of ~200,000 artificial stars have been simulated in each 
cluster in about 1000 simulation runs. Indeed, in order to not alter the crowding conditions, 
only 100-200 stars have been simulated in each run. The fraction of recovered objects in 
each magnitude interval was estimated as A = ■^ LS£ -, and a suitable completeness curve was 

■** sim 

obtained in each of the A,B,C regions (see right panels of Fig. 5). 

Star counts in each radial region have been finally corrected for incompleteness, by 
dividing each observed distribution by the corresponding A factors 2 . The total number of 
stars has been finally obtained by summing the completeness-corrected number of stars of 
the A,B,C regions. 

As already discussed in previous papers, the artificial star technique provides only a 
first-order correction. In fact, the observed distribution is, in principle, distorted because of 
two main phenomena: the loss of faint stars due to incompleteness and an excess of bright 
stars due to possible blending effects of two or more faint stars into a brighter one. Only 
the first effect is taken into account by the artificial star simulation. Blending from faint MS 
stars is a more complicated effect to simulate, however in the near IR it is negligible (see 



2 Note that the number of stars lost for incompleteness (n comp ) in each bin of magnitude is 
where n t, s is the number of stars observed in that bin. 
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also discussion in Testa et al. 1999). 

Another important effect which needs to be investigated, is the degree of contamination 
of the selected samples by the foreground/background stars. In this paper we have applied 
a statistical decontamination technique, using a control field adjacent to the clusters. The 
total number of stars observed in each evolutionary sequence (AGB, RGB and He-clump) 
has been counted accordingly to the selection boxes both in the cluster (see Figs. 1 and 3) 
and field (see Figs. 2 and 4) CMDs, and corrected for incompleteness (see above). The star 
counts in the field population have been scaled to take into account the different surveyed 
area, and their contribution have been subtracted from the cluster population. 

In summary, for each radial region, each selection box corresponding to each evolutionary 
stage has been divided in bins of magnitude (typically 0.2 mag wide). Then, the "corrected" 
number of stars in each bin has been computed as follows: 

where n b. s is the number of stars observed in that bin, n comp is the number of stars lost for 
incompleteness, n/ is the expected number of field stars. 

Analogously, the total luminosity of each evolutionary stage can be computed accord- 
ingly to the following relation: 



a 



comp 



-l^corr — \ / j L'i ) "•" \ n comp X i-i eq ) [Jlf X -L e q) 
i=l 

where the term Y17=i ^ bs * s ^ ne total luminosity of the stars observed in a given bin, n, 
is the number of stars lost for incompleteness, rif is the expected number of field stars, and 
L eq is the equivalent luminosity of that bin, that is the luminosity of a star with magnitude 
equal to the mean value of the bin. 

Finally, star counts and total luminosity of each evolutionary stage have been obtained 
by summing the contribution of all the bins in each selection box. 



5.2. Integrated magnitudes 

In order to properly perform cluster to cluster comparisons, one needs to take into 
account the size of the total cluster population. Hence both star counts and luminosities 
needs to be normalized to a reference population or to the cluster integrated luminosity. 
In previous papers we followed both the approaches (in F95 and F04 the RGB populations 
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were normalized to the cluster integrated luminosity taken from the literature, and to the 
He-clump population, respectively). 

The large FoV of SOFI offer the opportunity to independently determine the near- 
infrared integrated magnitudes for the program clusters. In doing this, we adopted a simple 
approach, by performing aperture photometry over the entire cluster extension (typically 
90"). In order to correct for the field contamination, an equivalent aperture photometry has 
been also performed on each control field and the resulting luminosity has been subtracted 
from the cluster value. 

The cluster center has been computed by applying a standard technique - Calzetti et 
al. (see for example 1993) - which uses the knowledge of the position of individual stars in 
the innermost region of the cluster, allowing a high precision determination of the center of 
gravity. Hence, by applying the procedure described in Montegriffo et al. (1995) we computed 
Cgrav by simply averaging the a and 5 coordinates of stars lying within a fixed radius from 
a first-guess center estimated by eye. The barycenter of the stars is then derived iteratively 
(see also Ferraro et al. (2003b)). The center of gravity (C grav ) of the programme clusters are 
listed in Table 1. Our new estimates turn out to be reasonably consistent (within ~ 10") 
with available determinations (as those in the SIMBAD astronomical database by the CDS, 
Strasbourg). The typical la uncertainty of our estimates is ~ 5 pixels corresponding to 
1.5'/ in both aj2ooo and <5j2ooo- The position of the center of 2 clusters, namely NGC2136 
and NGC2173, appear to be significantly (up to 2 arcmin) different from the SIMBAD 
coordinates. 

The case of NGC 2136 deserves an additional comment. Indeed, a small twin cluster, 
namely NGC 2137, is present at an angular distance of 1.34' (Hilker et al. 1995). Since its in- 
tegrated luminosity, although significantly fainter, is contaminating the aperture photometry 
of NGC 2136, it has been properly subtracted. 

Integrated K magnitude, colors and derived bolometric luminosities in the K band and 
in bolometric for the entire sample of 19 clusters are listed in Table 2. 



5.3. Error budget 

Formal errors are directly estimated from the photometric samples, by assuming that 
star counts follow the Poisson statistics. The errorbars for the various population ratios (by 
number and/or by luminosity) have been computed accordingly to the following formula 

y/R* -ai + vf, 

OR- ^ 
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with R = N/D, N being the numerator and D the denominator of the ratio. 

In the computation of the population ratios different error sources are at work, depend- 
ing on the observable. 

• total cluster luminosity: the main source of uncertainty in this case is the positioning 
of the cluster center. We estimate that an off-centering of 5 pixels corresponds to a 5% 
variation in luminosity. An additional uncertainty of ~ 10% has been considered in the 
computation of bolometric luminosities, in order to take into account the uncertainty 
in the bolometric corrections. 

• AGB luminosity: a conservative uncertainty of 0.2 mag in setting the faint end of the 
AGB luminosity distribution, implies a ~5% variation in the total AGB luminosity. 
However, for this observable the major source of uncertainty is the random error as- 
sociated to the number of detected AGB stars (in the Poisson regime a oc VNagb), 
which can suffer large fluctuations due to the small number statistics. On average, the 

tK 

overall cr R associated to the T jf B ratio turns out to be 30%. 

L TOT 

• number and luminosity of C-stars: as in the case of AGB stars, these observables and 
their associated errors suffer large fluctuations due to the small number statistics. On 

L K 

average, the overall Cr associated to the ^ star ratio turns out to be ~50%. 

• RGB luminosity: in the older clusters (see Fig. 3) the RGB is well populated, hence the 
estimated luminosity is much less affected by statistical fluctuations or by the selection 
box definition. On average, we estimate a a R ~20% for the RGB population ratios. 



6. The AGB and C-rich star contribution to the cluster light 

Theoretical models (Renzini & Buzzoni 1986; Maraston 1998, 2005) predict that the 
most important contributors to the integrated cluster light between 10 8 and 10 9 yrs are 
AGB stars. The AGB population includes both O-rich (M-type) and C-rich stars. During 
the thermal pulsing phase (hereafter TP-AGB) an AGB becomes C-rich if it undergoes the 
third dredge-up mixing process (see e.g. Iben & Renzini 1983). The presence of C-stars in 
stellar clusters depends on their age and metallicity (Renzini & Voli 1981). 

In intermediate age clusters the bulk of the AGB population is more luminous than 
the RGB Tip, and a minor overlap does exist between the faintest end of the AGB and 
the brightest portion of the RGB. Here we use our data set in order to investigate the 
contribution to the cluster luminosity of the brightest portion of the AGB populations as a 
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function of the cluster age. In order to consider the entire database that spans a large range 
of ages (from 80 Myr to 3.5 Gyr, see Table 1), in the following we consider only AGB stars 
brighter than K«i 12.3, which represents the RGB Tip level for the oldest clusters in our 
sample (see F04). 

The left panel of Fig. 6 shows the brightest portion of the Ko, (J — K)o cumulative 
CMD, where all the stars detected in the 19 surveyed clusters are plotted. The selection 
box adopted to sample the bright AGB population is over-plotted to the diagram. The right 
panel of Fig. 6 shows the cumulative (J — H) , (H — K) color-color diagram for the selected 
AGB stars. This diagram is especially suitable to isolate C-stars, since they are significantly 
redder than O-rich stars (see also Cioni et al. 2005). As shown in Fig. 6 a population of 26 
candidate C-stars (plotted as filled circles) has been identified on the basis of their extremely 
red colors in the 19 surveyed clusters. 

Note also that the artificial star experiments demonstrated that stars lying in the bright- 
est portion of the RGB can be safely recovered (with an overall photometric uncertainty of 
~ 0.03 mag) even in the innermost region of the clusters, excluding the possibility that 
blending of RGB stars could produce spurious bright objects lying within the AGB selection 
box. 

We used both the cumulative color- magnitude and color-color diagrams in Fig. 6 to make 
a census (both by number and luminosity) of the AGB stars brighter than the RGB Tip as 
well as of C-stars in each cluster. Although the number and the luminosity of AGB stars 
are subjected to large fluctuations due to the small-number statistics, we still performed a 
statistical decontamination, following the procedure described in Sect. 5.1. 

The number of AGB stars counted in each clusters and the number adopted after the 
field de-contamination are listed in Table 3. Once the accurate census of the AGB population 
(by number and luminosity) is available for all the sampled clusters, a number of suitable 
diagnostics tools can be used in order to study the AGB properties as a function of the age. 

The upper panel of Fig. 7 shows the ratio between the AGB and the cluster integrated 
K-band luminosity, as a function of the cluster age. It is remarkable the rapid increase (up 
to a factor 2) of the AGB luminosity at rs 200 Myr which reaches its maximum contribution 
in the 300-700 Myr range followed by a rapid decrease. Note that in 2 clusters (namely 
NGC2108 and NGC1987 at s = 35 - 36 corresponding to t ~ 600 - 700Myr) the brightest 
portion of the AGB population account for m 90% of the total cluster luminosity. These 
results are in good agreement with F95, who found that the maximum contribution of the 
AGB to the cluster light occurs at s = 35, corresponding to an age of ~ 600 Myr. 

The lower panel of Fig. 7 shows the same ratio as in the upper panel but with the 
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clusters grouped in five age bins accordingly to their s-parameter, namely s=23-26, 27-31, 
35-36, 37-39 and 40-45, respectively. For each bin we computed the weighted mean and 
the corresponding standard deviation. Theoretical predictions from Maraston (1998, 2005) 
for [Z/H]=-0.33 are also plotted. In these models the TP-AGB energetics was calibrated 
with previous (FM90, F95) intermediate age MC cluster data (see Maraston (1998) for full 
details). Our new observations nicely confirm the modelling and those early results. Old 
canonical models of stellar evolution (Renzini & Buzzoni 1986) were dating the occurrence 
of the AGB Ph-T by Renzini & Buzzoni (1986) at significantly earlier epochs (~ 10 7 yrs) 
of the stellar lifetime with respect to the new models. This discrepancy is due to a different 
treatment for the TP-AGB stars that experiencing the envelope burning process (Renzini & 
Voli 1981; Bloecker & Schoenberner 1991), as widely discussed in Maraston (1998). 

As a further evidence, Fig. 8 shows the ratio between the number of C-stars and total 
bolometric luminosity and the ratio between the K luminosity of C-stars and the total K- 
band cluster luminosity, as a function of the cluster age. The number of C-stars in each 
cluster is listed in Table 3. We assumed all these stars being cluster members. Indeed, we 
estimate that the probability to find a field C-star within the sampled cluster area is <30%. 
The contribution of C-stars to the total cluster luminosity as a function of the cluster age, 
closely follows the one shown by the entire AGB population (see Fig. 7) and it turns out 
to be larger than 50% in the 700 - 1000 Myr age range. In NGC 2190 (with s = 36, 
hence t = 730 Myr) the C-stars accounts for 70% of its total luminosity. Previous works by 
FMB90 also found that the fraction of luminosity from bright AGB and C-stars is maximum 
for SWB V clusters (i.e. s between 35 and 40). 

It is worth noticing the case of NGC 2249, whose age (log t = 8.72) corresponds to 
the epoch when the AGB contribution is expected to reach its maximum. Conversely, both 
Figs. 7 and 8 shows that NGG 2249 (marked with an open circle) has a very low AGB 
luminosity for its age. Indeed no C-stars and only 1 AGB have been detected in this cluster. 
On the other hand, NGC2249 is the least luminous cluster in our sample (Lk ~ 5 x 10 4 L q ), 
hence the fastest evolutionary stages (as the AGB) are expected to be intrinsecally poorly 
populated in its CMD. 



7. The RGB Ph-T 

For the analysis of the RGB Ph-T we have considered the clusters in our data-base 
with s > 34 (the 6 oldest clusters in the sample presented here and the 6 discussed in F04). 
Obviously the 7 youngest clusters in our sample (Fig. 1) are not considered in the following 
discussion because they have not developed a populous RGB yet. The considered sample 



-14- 



covers a wide range of ages (from 500 Myr up to 3.5 Gyr) significantly extending the age 
range covered by F04 and it allows to probe the entire developing of the RGB. 

In order to compute the RGB population ratios we have adopted the same procedure 
described in F04. Using the cumulative K Q , (J — K) CMD for the 12 clusters we identified 
the mean loci of the upper RGB and the He-Clump and defined the corresponding selection 
boxes sampling these populations (see F04). As discussed in Sect. 6 the RGB Tip is expected 
to be at K Q ^12.3. As an example, Fig. 9 shows the de-reddened CMD of NGC 1651 where 
the two selection boxes for the RGB and He-clump population, respectively, have been 
plotted. 

Since the photometric errors can significantly broaden the faint sequences, the size of 
the boxes including the base of the RGB and the He-clump has been conservatively assumed 
to be ~ 5 times the photometric error at that level of magnitude 3 . 

Population counts and luminosities for stars in the RGB and He-clump evolutionary 
stage have been obtained and corrected for incompleteness and field contamination accord- 
ingly with the procedure discussed in Sect. 5.1. The results are listed in Table 4 and plotted 
in Figs. 10 and 11. Fig. 10 shows the behavior of the number of RGB stars normalized to 
the number of He-Clump stars as a function of the cluster age. Fig. 11 shows the bolometric 
luminosities of RGB stars normalized to the bolometric luminosities of He-Clump stars as 
a function of the cluster age (upper panel) and the bolometric luminosities of RGB stars 
normalized to the bolometric luminosities of entire cluster as a function of age (lower panel). 

Note that, in a few high density clusters (see Table 4) with severe crowding, completeness 
drops down to 60% at the He-clump magnitude level in the innermost A region (see Fig. 5). 
Hence, in these clusters star counts and luminosities have been computed only in the outer 
B and C regions. 

As already discussed in F04, at an age of ~ 500 Myr the rapid increase of the RGB 
population ratios (by a factor of ~3 in number and ~4 in luminosity) in a timescale as short 
as m 400 Myr flags the occurring of the RGB Ph-T. At the age of ~ 900 Myr a progressive 
flattening of the ratios suggests that the full developing of an extended and well populated 
RGB has occurred. The overall increase of the population ratios between m 500 Myr and 
~ 3.5 Gyr is a factor ^5 by number and ~7 by luminosity. 

Empirical data have been compared to theoretical predictions. In Fig. 10 and Fig. 11, we 
over-plotted the predictions of canonical models (Maraston 1998, 2005) with [Z/H] = —0.33. 



3 Note that, since the boxes sample the bulk of the population along each evolutionary stage, a slightly 
different assumption in the selection box size has a negligible impact on the overall results. 
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The models nicely agree with the observations over the entire range of considered ages, well 
describing the epoch, the duration and the increasing contribution of the RGB phase. 

It is worth noticing the population ratio excess (both by number and by luminosity, see 
Figs. 10 and 11) of NGC 1783, when compared with other clusters with similar values of the 
s— parameter (i.e., NGC 2231). It is likely that this cluster be older than suggested by the 
s— parameter, since its CMD (see Fig. 3) shows a fully populated RGB, typical of clusters 
with s-parameter>40. Indeed, its RGB morphology is more similar to that one of clusters 
such as NGC 1806, NGC 2173 and NGC 1978 rather than that one of NGC 2231 (see Fig. 2 
in F04). This evidence further supports the urgency of a new homogeneous calibration of 
the age scale of LMC clusters. 



8. Conclusions 

We have used our extensive near IR database of 19 LMC clusters spanning a wide range 
of ages (from 80 Myr to 3.5 Gyr) in order to investigate the developing in time of the AGB 
and RGB evolutionary stages. 

The behavior of the quantitative contribution to the total cluster luminosity of both 
the RGB and AGB in terms of age has been investigated and turns out to be nicely in 
agreement with theoretical prediction of canonical models. The AGB contribution to the 
total cluster luminosity starts to be significant at ~ 200 Myr and reaches its maximum peak 
at 5 — 600 Myr (see Fig. 7), when the RGB Ph-T is starting (see F04 and Figs. 10 and 11). 
Both events are sharp and last a few hundreds Myr only. 

The epoch and duration of the AGB phase-transition estimated from our data con- 
firm the semi-empirical modelling of Maraston (1998, 2005). The latter was based on the 
calibration of the AGB-fuel consumption of canonical tracks, and its partition between C 
and O-type stars, with previous MC GCs data (FMB90, F95). Furthermore, the epoch and 
duration of the RGB phase transition derived from our data agree very well with the the- 
oretical predictions of simple stellar population models based on canonical tracks and the 
fuel-consumption approach (Maraston 2005). 

This study represents a further step towards the calibration of the integrated properties 
of SSPs as a powerful diagnostics for dating complex stellar populations like galaxies. 

The financial support by the Ministero dell'Istruzione, Universita e Ricerca (MIUR) is 
kindly acknowledged. 
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Table 1. Main parameters of the entire sample of observed LMC clusters. 



Cluster 


a (J2000) 


5 (J2000) 


s 


age (Myr) 


[Fe/H] 


E(B-V) 




NGC 2164 


05:58:55.65 


-68:31:00.75 


23 


81 


-0.60" 


0.10 


this paper 


NGC 2157 


05:57:36.74 


-69:11:53.58 


25 


114 


-0.60" 


0.10 


this paper 


NGC 2136 


05:52:58.54 


-69:29:32.32 


26 


135 


-0.55 6 


0.10 


this paper 


NGC 2031 


05:33:39.00 


-70:59:14.54 


27 


160 


-0.52 6 


0.18 


this paper 


NGC 1866 


05:13:38.88 


-65:27:53.30 


27 


160 


-0.50 d 


0.10 


this paper 


NGC 2134 


05:51:57.54 


-71:05:51.63 


28 


190 


-1.00" 


0.10 


this paper 


NGC 1831 


05:06:16.47 


-64:55:12.76 


31 


315 


+0.01 c 


0.10 


this paper 


NGC 2249 


06:25:49.50 


-68:55:14.25 


34 


524 


-0.12" 


0.10 


F04 


NGC 1987 


05:27:17.29 


-70:43:56.78 


35 


620 


-1.00" 


0.12 


F04 


NGC 2209 


06:08:34.87 


-73:50:06.46 


35 


620 


-1.20" 


0.07 


F04 


NGC 2108 


05:43:57.30 


-69:10:55.93 


36 


734 


-1.20" 


0.18 


F04 


NGC 2190 


06:01:00.67 


-74:43:29.10 


36 


734 


-0.12 c 


0.10 


F04 


NGC 2231 


06:20:43.67 


-67:31:13.05 


37 


869 


-0.67 c 


0.08 


F04 


NGC 1783 


04:59:08.42 


-65:59:12.75 


37 


869 


-0.45" 


0.10 


this paper 


NGC 1651 


04:37:33.86 


-70:35:09.24 


39 


1218 


-0.37 c 


0.10 


this paper 


NGC 2162 


06:00:30.20 


-63:43:15.27 


39 


1218 


-0.23 c 


0.07 


this paper 


NGC 1806 


05:02:11.87 


-67:59:10.11 


40 


1442 


-0.23 c 


0.12 


this paper 


NGC 2173 


05:57:59.28 


-72:58:42.83 


42 


2021 


-0.24 c 


0.07 


this paper 


NGC 1978 


05:28:45.34 


-66:14:09.12 


45 


3353 


-0.96 d 


0.10 


this paper 



Note. 



The s-parameter is from Elson & Fall (1985) and Girardi et al. (1995). Metallicity are 
from (a) Sagar & Pandey (1989), (b) Dirsch et al. (2000), (c) Olszewski et al. (1991), and 
(d) Hill et al. (2000). Reddening is from Persson et al. (1983). 
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Table 2. Integrated K-magnitude, colors and luminosities of the target clusters. 



Cluster 


s 


K 


J-K 


H-K 


jK 
'-'TOT 


j Bol 

±J TOT 


NGC 2164 


23 


8.575 


0.606 


0.122 


19.76 


12.82 


NGC 2157 


25 


8.286 


0.658 


0.195 


25.79 


15.18 


NGC 2136 


26 


7.968 


0.674 


0.125 


34.57 


19.78 


NGC 2031 


27 


8.219 


0.936 


0.323 


28.14 


11.36 


NGC 1866 


27 


7.282 


0.743 


0.155 


65.02 


33.13 


NGC 2134 


28 


9.161 


0.742 


0.187 


11.52 


5.88 


NGC 1831 


31 


8.313 


0.802 


0.231 


25.16 


11.67 


NGC 2249 


34 


9.983 


0.934 


0.342 


5.40 


2.06 


NGC 1987 


35 


8.811 


0.988 


0.316 


16.00 


5.75 


NGC 2209 


35 


8.955 


1.224 


0.345 


13.79 


3.54 


NGC 2108 


36 


8.807 


1.148 


0.403 


16.38 


4.96 


NGC 2190 


36 


9.269 


1.193 


0.321 


10.43 


2.83 


NGC 2231 


37 


9.329 


1.096 


0.328 


9.81 


2.97 


NGC 1783 


37 


7.091 


1.035 


0.282 


77.53 


25.77 


NGC 1651 


39 


8.895 


0.973 


0.301 


14.72 


5.32 


NGC 2162 


39 


9.071 


1.253 


0.372 


12.40 


3.07 


NGC 1806 


40 


7.076 


1.055 


0.271 


79.11 


25.97 


NGC 2173 


42 


9.05 


1.033 


0.297 


12.33 


4.02 


NGC 1978 


45 


7.185 


0.894 


0.288 


71.10 


28.69 



Note. - 
Integrated magnitudes are decontaminated for the field 
contribution. K-band and bolometric total luminosities 
are in units of 1O 4 L . 
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Table 3. Star counts and luminosities for AGB and C-stars 



Cluster 


s 


N^b 


]\jField 
ly AGB 


JVfDec 
ly AGB 


T K 

U AGB 

T K 

^TOT 


Nc-star 


NC-star 
J Bol 
^TOT 


Nc-star 

T K 


T K 

^C-star 

J K 

^TOT 


NGC 2164 


23 


2 





2 


0.10 














NGC 2157 


25 


9 


1 


8 


0.52 














NGC 2136 


26 


9 


1 


8 


0.31 


1 


0.05 


0.03 


0.09 


NGC 2031 


27 


7 


1 


6 


0.64 














NGC 1866 


27 


12 





12 


0.36 














NGC 2134 


28 


2 


1 


1 


0.71 














NGC 1831 


31 


7 


1 


6 


0.45 


3 


0.26 


0.12 


0.31 


NGC 2249 


34 


1 





1 


0.32 














NGC 1987 


35 


9 


4 


5 


0.90 


3 


0.52 


0.19 


0.61 


NGC 2209 


35 


4 





4 


0.72 


2 


0.56 


0.14 


0.58 


NGC 2108 


36 


5 


1 


4 


0.87 


1 


0.20 


0.06 


0.34 


NGC 2190 


36 


2 





2 


0.74 


2 


0.71 


0.19 


0.72 


NGC 2231 


37 


1 





1 


0.32 


1 


0.34 


0.10 


0.32 


NGC 1783 


37 


16 


1 


15 


0.36 


2 


0.08 


0.03 


0.11 


NGC 1651 


39 


4 





4 


0.51 


1 


0.19 


0.07 


0.09 


NGC 2162 


39 


4 


1 


3 


0.73 


1 


0.32 


0.08 


0.59 


NGC 1806 


40 


13 


4 


9 


0.22 


4 


0.15 


0.05 


0.17 


NGC 2173 


42 


5 


1 


4 


0.53 


1 


0.25 


0.08 


0.15 


NGC 1978 


45 


13 


1 


12 


0.25 


4 


0.14 


0.06 


0.13 



Note. - 
Star counts are corrected for incompleteness. K-band and bolometric luminosities are 
in units of 10 4 L . 
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Table 4. Star counts and bolometric luminosities for RGB and He-Clump stars. 



Cluster 


s 


N RGB (a) 


N He -c (a) 


T bol (b) 
L RGB 


T bol (b) 
L He-C 




NGC 2249 


34 


9 


98 


0.16 


0.51 


F04 


NGC 1987 


35 


42 


322 


0.92 


2.05 


F04 


NGC 2209 


35 


24 


160 


0.61 


0.81 


F04 


NGC 2108 


36 


40 


231 


1.11 


1.38 


F04 


NGC 2190 


36 


28 


174 


0.94 


0.89 


F04 


NGC 2231 


37 


36 


114 


0.71 


0.59 


F04 


NGC 1783^ 


37 


150 


352 


4.58 


1.98 


this paper 


NGC 1651 


39 


43 


177 


1.16 


1.09 


this paper 


NGC 2162 


39 


40 


143 


0.99 


0.78 


this paper 


NGC 1806 (c) 


40 


75 


218 


2.25 


1.30 


this paper 


NGC 2173 


42 


36 


84 


1.11 


0.48 


this paper 


NGC 1978 (c) 


45 


182 


402 


5.03 


2.23 


this paper 



Note. - 

(a) Star counts are corrected for incompleteness and field contamination. 

(b) Bolometric luminosities are in units of 10 4 L . 

(c) Due to severe crowding conditions, star counts and luminosities have 
been computed only in the outer B and C regions, see Fig.5. 
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Fig. 1.— Observed (K, J-K) CMDs of the 7 observed LMC clusters with s = 23-31. In the 
last panel a sketch of the CMD loci dominated by the cluster (grey regions) and the LMC 
field (dashed box) populations are also shown for sake of clarity. 
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Fig. 2. — Observed (K, J-K) CMDs of the fields adjacent to the 7 observed LMC clusters 
with s = 23 — 31. 
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Fig. 3. — Observed (K, J-K) CMDs of the 6 observed LMC clusters with s = 39 — 45. 
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Fig. 4. — Observed (K, J-K) CMDs of the fields adjacent to the 6 observed LMC clusters 
with s = 39 - 45 
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Fig. 5. — An example of cluster radial mapping. Left panel: the cluster frame is divided into 
3 concentric annuli, to account for different crowding conditions. Right panel: completeness 
curve for each radial sub-region, as labeled in the left panel. 
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Fig. 6. — Le/t panel: cumulative, de-reddened K , (J — K) CMD for the entire cluster 
sample. The selection box adopted to isolate the AGB population (large circles) is shown. 
Right panel: de-reddened color-color (J — H) , (H — K) diagram of the AGB stars. In both 
panels open circles are O-rich AGB, filled circles are C-stars. Solid line boxes to distinguish 
C-stars and Long Period Variables (LPV) are from Bessell & Brett. (1988), (see also Ferraro 
et al. 1995). The mean locus for K giant stars (solid line) is from Frogel et al. (1978). 
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Fig. 7. — Upper panel: observed AGB contribution to the total cluster K-band luminosity as 
a function of age. The open circle marks the intrinsically poor populated cluster NGC 2249. 
Lower panel: weighted mean and standard deviation of the same ratio with the clusters 
grouped in five age bins, accordingly with their s-parameter, namely s=23-26, 27-31, 35- 
36, 37-39 and 40-45, respectively (the cluster NGC 2249 has been excluded). Theoretical 
predictions for the temporal evolution of the entire AGB (E-AGB and TP-AGB, solid line) 
and for the dominant TP-AGB (dashed line) are overplotted. Both models are computed at 
[Z/H]= -0.33. 
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Fig. 8. — Upper panel: the number of C-stars normalized to the K-band luminosity of the 
cluster as a function of age. Lower panel: the K-band luminosity of the C-stars normalized 
to the total cluster luminosity as a function of age. The open circles marks the intrinsically 
poor populated cluster NGC 2249. 
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Fig. 9. — An example of a K , (J — K) de-reddened CMD with the selection boxes adopted 
to distinguish the RGB and the He-Clump populations for clusters with s > 39 (see Fig. 3). 
The position of the RGB Tip is also indicated. 



-31- 



0.5 - 



0.4 - 



u 



* 0.3 - 



O 
OS 



^ 



0.2 - 



0.1 - 




1000 



2000 

Age (Myr) 



3000 



4000 



Fig. 10. — Ratio between the number of the bright RGB and He-Clump stars as a function of 
age for the 12 clusters with s > 33: filled circles are the 6 clusters with s = 39 — 45 presented 
in this paper (see Fig. 3), filled squares are the 6 clusters from F04. Stars belonging to 
the two populations are selected accordingly to the selection boxes shown in Fig. 9. The 
solid line represents the prediction of the canonical theoretical model with [Z/H]= —0.33 
(Maraston 1998). 
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Fig. 11. — Upper panel: the bolometric luminosity of the RGB normalized to the He-Clump 
as a function of age for the 12 clusters with s > 33. Symbols and lines are as in Fig. 10. 
The line represents the prediction of the canonical theoretical model with [Z/H]= —0.33 
(Maraston 1998). Lower panel: the bolometric luminosity of the RGB normalized to the 
total bolometric luminosity for the 12 clusters. The line represents the prediction of the 
canonical theoretical model with [Z/H]= —0.33 (Maraston 2005). 



